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THE  EDGE-FUNCTION  METHOD  (E.F. M . ) 


ABSTRACT  OF  INVITED  LECTURE  by  P.  M.  QUINLAN  to 

The  International  Symposium  on  Innovative  Numerical  Analysis  in  Applied 
Engineering  Science,  Versailles,  23-27  May  1977. 

(The  Research  was  carried  out  in  cooperation  with  Professor  J.  E.  Fitzgerald 
and  Professor  S.  N.  Atluri,  both  of  Georgia  Institute  of  Technology). 
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SCOPE  OF  THE  EDGE-FUNCTION  METHOD 

The  Edge-Function  Method  as  applied  to  two-dimensional  problems  is  a pro-  i 
cedure  for  obtaining  exact  solutions  to  acceptable  Mathematical  Models  of  ; 

practical  boundary  value  problems  in  regular  or  irregular  polygonal  regions 
with,  or  without,  cavities.  The  method  was  originally  developed  by  Quinlan  [].)  [ 
for  the  torsion  of  prismatic  bars  of  polygonal  cross-section.  i 

In  the  succeeding  12  years  it  has  been  successfully  applied  to  problems  ’ 

in  the  bending  of  isotropic  thin  plates  [2,3,4];  coupled  linear  systems  in  j 

elastostatics  [4, 3, 6,7]  and  in  moderately  thick  plates  and  shallow  shells  [8]; 
cracks  and  stress  concentrations  in  elastostatics  [9];  vibrations  of  Thin 
Plates  [10]  and  vibrations  of  Shallow  Shells  [11].  Selections  from  relevant  i I 
abstracts  are: 

* ' ^ '*  J 

"A  general  computer  method  is  presented  of  solving  a very  wide  range  of  ' 
swept-back  wing  and  skew  bridge  problems'- with  any  boundary  conditions  and  load- 
ings, with  or  without  cut-outs  and  elastic  column  supports".  [2] 

"The  Edge-Function  Method  may  be  described  as  a piecing  together  of 
"asymptotic"  solutions  to  a set  of  linear  p.d.e.'s  for  the  several  parts  of  a i 
domain  D to  satisfy  the  boundary  conditions  in  a discrete  least  squares  sense.  | ! 

In  the  present  paper,  application  is  made  to  problems  in  elastostatics,  where  D i 
may  be  multiply  connected, 'bounded  externally  by  an  irregular  polygon  and  with  ! j 
curved  holes  in  the  interior.  The  various  "asymptotic"  problems  are  solved  for  . 

(1)  the  infinite  sector,  for  Vertex  Functions,  (2)  the  semi-infinite  region, 
leading  to  Edge-Functions,  and  (3)  a curved  hole  in  an  infinite  region,  leading  , 
to  Mapped  Polar  Biharmonics . [7] 

"A  new  method  for  the  numerical  solution  of  Keissner-type  plate  and 
displacoment-typc  snallow  shell  equations  is  developed.  By  deriving  for  an 
arbitrary  edge  of  a polygonal  plate  or  shallow  spl  rical  or  rectangular 
cylindrical,  hyperbolic  paraboloidal,  or  elliptic  nraboloidal  shell  a set  of 
"edge-functions",  it  is  possible  to  obtain  accurate  solutions  to  these  j 

structural  components  subjected  to. various  support  and  load  conditions".  [8]  j 
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"The  present  investigation  extends  the  method  o£  Edge-Functions  to  the 
determination  of  natural  frequencies  and  associated  mode  shapes  of  free 
vibration  of  thin  elastic  plates  with  a variety  of  boundary  conditions.  The 
Edge-Function  technique  essentially  associates  an  independent  coordinate  system 
with  every  edge  of  the  plate,  and  for  every  edge  employs  functions  (satisfying 
the  governing  plate  vibration  equations)  that  rapidly  decay  with  increasing 
distance  from  the  plate  boundary.  Appropriate  superposition  of  solutions 
stemming  from  all  edges  of  the  plate  leads  to  an  approximate  representation  of 
the  time-dependent  deflection  surface".  [10] 

"It  can  confidently  be  concluded  that  the  edge-function  method,  using 
Edge-Functions  and  the  newly  developed  "Shell  Polar"  functions,  provides  a 
rapid  means  of  determining  accurate  values  for  the  frequencies  of  free  vibration 
and  mode  shapes.  It  computes,  as  a routine,  approximate  root  mean  square 
boundary  residuals  for  each  boundary  function,  thus  providing  a practical 
measure  for  judging  the  acceptability  of  the  solution  proferred".  [11] 

FUNDAMENTALS  OF  E.F.H.  VIA  LAPLACE'S  EQUATION 

Section  2 of  the  paper  aims  to  illustrate  as  simply  as  possible  the  main 
algebraic  and  programming  features  of  by  applying  it  to  Laplace's 

equation,  V^u  - 0,  for  polygonal  regions  with  elliptical  cavities. 

On  superposing  1 harmonic  solutions,  u.,  that  represent  the  main 
characteristics  of  the  problem,  we  obtain 

u - Z A u.  ; V2u  = 0 , ~ (1) 

r=l  11  1 'S* 

where  the  I superposition  constants  A.-ar'e  evaluated  numerically  to  satisfy,  in 
an  approximate  way,  the  specified  boundary  conditions. 

The  method  of  harmonic  fitting  is  developed,  by  expanding  the  boundary 
identity,  say  u=0,  on  a typical  side  j,  0<x'<a’,  in  a Fourier  sine  series 

OO  1 . , 


0 H u(x') 


z c.. 

N=1  N 


C = — , j“  u(x')sin  nx'dx’ 
N a 0 


The  series  cau,  in  practice, /be'-  truncated  at  N'  terms  giving  the  approximation 

N ' I 

0 = u(x')  ~ I C sin  nx’  ; 0 < x'  < a'  ; n = Nr/a’  (3) 

N=1  " 

3 

provided  that  the  series  has  convergence  of  order  1/N  to  ensure  continuity  of 
u(x'),  which  requires  the  two  point,  or  vertex,  equations 

u(x')  *=  0 ; x'  = 0 and  x*  = a'  (4) 

Approximation  (3)  then  leads  to  the  harmonic  equations 

CN  = 0 ; N = 1 N’  (5) 

On  writing  a pair  of  vertex  equations  and  N'  harmonic  equations  for  each  side  ofj 
the  polygon,  a system  of  I linear  equations  results  for  the  unknowns  A. , which 
equations  may  be  written  as 

K A = b , _ (6) 

where  K is  the  coefficients  matrix,  and  b the  right  hand  column  vector. 

The  integral  for  in  (2)  may  require  numerical  integration.  A better 
approach  is  to  interpret  (3)  as  the  fitting  of  a trigonometric  expression  to 
u(x'),  and  on  using  a discrete  least  squares  criterion  of  fit,  we  obtain 

cn  ’ TT7T  u(xM)sin  nxA  ; M’  - N'  (7) 

where  x^  are  M’  interior  equidistant  points  on  side  x'.  Note  that  the  series  J 
for  is  formally  obtained  by  evaluating  integral  (3)  numerically  by  the 
trapezoidal  rule  with  M'+l  segments. 

The  basic  problems  - for  the  hasic  regions:  infinite  sector,  somi- 
infitiitc  strip  and  a hole  in  an  infinite  plane  - encountered  in  putting  together! 
a suitable  solution  mix  (1),  are  then  solved.  The  resulting  functions,  called 
basic  functions,  arc  called  the  Vertex,  Edge  and  Mapped  Polar  functions  \ 
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i 
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respectively.  Extension  is  made  to  mixed  boundary  conditions,  and  to  polygons 
with  elliptical  indentations  by  the  use  of  conformal  mapping.  A simple  computer 
program  is  available  to  illustrate  the  basic  features  of  E.F.M.  for  V^u  = 0 for 
a polygonal  region,  when  u is  specified  on  all  boundary  segments. 

Fourier  Cosine,  and  full  Fourier  series,  arc  developed  as  alternatives  to  j 
expansion  (2),  as  these  are  required  when  dealing  with  Elastostatics  and 
Orthotropic  Thin  Plates.  Finally  extension  is  made  to  a region  divided  into  a 
series  of  elements,  as  in  Finite  Element  work,  and  to  cases  where  singularities  : 
occur  either  within  or  on  the  boundaries  of  the  region.  ' 

A comprehensive  computer  program  for  Laplace/Foisson  problems  is  available 
with  illustrative  examples,  to  act  as  a program  prototype  for  the  E.F.M.  when 
applied  to  other  areas  - elastostatics,  thin  plates,  vibrations  etc.  i 

Recent  Research  A largely  self-contained  paper  [9]  is  due  for  publication.  This! 
sets  the  basis  for  dealing  with  coupled  linear  systems  through  the  medium  of 
elastostatics,  and  follows  the  general  lines  of  section  1 to  obtain  the  basic  ; 

functions.  Examples  are  given  involving  fracture  and  reentrant  angles  and  a 1 

comprehensive  program  is  available. 

A companion  work  [14]  on  isotropic  plates  updates  and  extends  considerably' 
the  earlier  work  [2],  while  [15]  a Ph.D.  thesis,  to  be  represented  in  June,  does' 
a comprehensive  job  on  aelotropic  and  orthotropic  thin  plates  for  polygonal 
regions  with  or  without  elliptic  holes.  This  latter  presents  a general  method 
for  determining  vertex  functions  directly  from  their  eigenvalue  determinant.  : 

I 

THREE-DIMENSIONAL  ELASTICITY  ' j 

Section  2 is  a progress  report  on  the  development  of  E.F.M.  for  three-  : 
dimensional  isotropic  elasticity.  The  Papkovich-Keuber  potentials  [15]  for  the.  • 
class  of  displacements  _ __ 

s = B + VX  ; V B = 0,.  (8)  ! 

are  the  key  to  obtaining  (1)  Plane  Functions  that  decay  exponentially  inwards 
from  each  plane  boundary  area,  and  (2)  Kcdge  Functions  that  model  the  singular 
behaviour  between  two  intersecting  plane  faces  - the  respective  analoges  of  the  ' 
Edge  and  Vertex  Functions  in  [9].  1 

The  boundary  conditions  on  each  plane  face  (x',y')  lead,  analogous  to  (2),  ; 
to  a plane  boundary  identity  i 

V (x\y’>  = 0,  , (9)  ' 

and  a discrete  lease  squares  surface  fitting  criterion  gives  a double  series  for 
harmoni'  fitting,  an-ilogoes  to  (7).  Thus  the  functions  and  methods  for  2-D  have 
their  counterparts  in  3-D  and  comparative  results  are  expected  before  the 
Symposium  for  the  problems  solved  by  Deak  [16]  using  Boundary  Integral. 

The  coef  f icic..ts  matrix  K arising  in  (6)  can  be  decomposed  into 

V - b - K2A,  (10) 

where,  since  the  terms  K,  are  of  minor  significance,  their  effects  can  be 
obtained  by  iteration  on' the  banded  matrix  for  kth-iteration: 

KjA  = b(k)  (11) 

Usually  three  iterations  provide  an  acceptable  solution,  the  time  spent  on  the 
iterations  being  only  a few  percent  of  that  required  to  solve  the  banded  matrix 
K^.  Thus  E.F.M.  is  now  firmly  within  the  banded  matrix  field,  and  is  capable  of 
dealing  with  any  number  of  connected  elements  as  in  Finite  Element  work. 

MAIN  ADVANTAGES 

The  effective  algorithmic  method  that  has  been  developed  for  controlling 
the  algebra  [9],  allied  to  a systems  approach  adopted  to  programming, 
facilitates  the  extension  of  E.F.M.  into  new  areas. 

Solutions  for  up  to  four  different  sizes  of  matrix  arc  computed  as  <1 
routine.  Each  solution  vector  provides  an  exact  solution  to  the  specified 
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problem  but  with  slightly  different  boundary  values,  which  are  computed  as  a 
routine,  the  differences  being  termed  the  boundary  res idnals . If  the  residuals 
are  within  the  limits  within  which  an  engineer  can  specify  the  boundary  values, 
then  E.F.M.  offers  as  "exact"  a solution  as  can  be  obtained  to  the  problem. 

Derived  quantities  e.g.  moments,  shears,  are  obtained  as  accurately  as  the 
nrimarv  quantity,  and  are  quickly  computed  for  several  truncation  levels. 

Perhaps  the  greatest  advantage  of  E.F.M.  is  that  no  pre-computer  process- 
ing of  the  problem  is  required.  Only  the  geometrical,  load  and  material  data  are 
required  and  an  E.F.M.  program  can  be  operated  as  a "black  box"  without  tha  user 
having  any  knowledge  of  E.F.M. 

ACKNOWLEDGEMENTS 

P.  M.  Quinlan  gratefully  acknowledges  research  grants  from  E.O.A.R.  from  1957-70, 
and  travel  grants  1972-76,  without  which  E.F.M.  would  never  have  been  developed. 

All  three  authors  are  pleased  to  acknowledge  a recent  grant  by  N.S.F.  Eng. 76  16413 
to  Georgia  Institute  of  Technology  for  research  on  "The  Edge-Function  Method  for 
3-Dimensional  Elasticity". 

REFERENCES 

1.  Quinlan,  P.  M.  (1964),  "The  Torsion  of  an  Irregular  Polygon",  Proc.Roy.Soc. , 

Vo 1 . 2S2A. 

2.  Quinlan,  P.  M.  (1968),  "Polygon  and  Swept-Back  plates  with  Cut-Outs  and 
Column  Supports",  Proc.OAR  Research  Applications  Conference,  Washington  1968. 

3.  Quinlan,  P.  M.  (1965),  "The  X Method  for  Skew  Plates",  Proc.  Fourth.  U.S. 

Mat.  Congr.  Ann.  Mech. 

4.  Quinlan,  P.  M.  (1971)  , "Final  Scientific  Report  on  Grant  AF  EOAR-69-0049", 
With  Computer  Program  QUINP,  as  an  appendix. 

5.  O'Callaghan,  M.  J.  (1972),  "The  Edge^Tuncticn  Method  for  Linear  Boundary 


Value  Problems  with  Discontinu 


Ph.B.  Thesis,  U.C.C. 


Quinlan,  P.  M.  (1974),  "The  Edge-Tune cion  Method  in  Elastostatics"  paper  in 
"Studies  in  Numerical  Analysis" ,, 'Academe  Pres:;.,  London  1974. 

Quinlan,  P.  M.  (1974),  "The  Edgcp^F unction  Method",  16th  British  Colloquium 
Theoretical  Mechanics,  St.  Andrews. 

Desmukh,  R.  S.  (1973),  "The Qylge-Function  Method  Applied  to  Moderately 
• Thick  Plates  and  Shallow  Shtflls",  Ph.D.  thesis.  University  of  Massachusetts. 
Quinlan,  P.  M.  and  O'Callaghan,  M.  J.  (1976),  "The  Edge-Function  Method  for 
Cracks  and  Stress  Concentrations  in  2-0  Problems",  to  be  published  in  Int.J. 
Nnm,  Math.'  Eng. 

Tai,  I.  H.  and  Nash,  W.  A.  (19/3),  "Vibrations  of  Thin  Plates  - A New 
Approach",  AFOSR  TR-74-0789,  University  of  Massachusetts. 

O'Callaghan,  M.  J. , Nash,  W.  4.  and  Quinlan,  P.  M.  (1975),  "Vibration 
Analysis  of  Thin  Shallow  Spherical  Shells  Using  the  Edge-Function  Method", 
AFOSR  TR-72-2340,  University  of  Massachusetts. 

Quinlan,  P.  M.  (1965),  "The  X Method  for  Rectangular  Plates",  Proc.Roy.Soc. , 
Vol . 288A. 

Quinlan,  P.  M.  and  Cusens,  A.  J.  (1977),  "The  Edge-Function  Method  in 

Bridge  Slab  Analysis",  (in  preparation). 

Grannell,  J.  (1977),  "The  Edge-Function  If  thod  for  Thin  Aeolotropic  Plates", 
Ph.  D.  Thesis,  U.  C. , Cork,  (to  bn  submitted  June  1977). 

Lurie,  A.  I.,  "Three-Dimensional  Problems  of  the  Theory  of  Elasticity", 

Wiley,  (1964). 

Deak,  Alexander,  L.  (1972),  "Numerical  Solution  of  3-D  Elasticity  Problems 
for  Solid  Rocket  Grains",  AFSC  Report  No.  AFRPL-TR-7 1-140 , Vol.  1. 

NOTE:  The  complete  paper  will  be  published  in  The  Internation- 

al Journal  of  Solids  and  Structures. 


